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ABSTRACT 

We report on X-ray observations made with the Rossi X-ray Timing Explorer of the black hole candidate (BHC) 
transient XTE J 1650-500 at the end of its first, and currently only, outburst. By monitoring the source at low 
luminosities over several months, we found 6 bright ^100 s X-ray flares and long time scale oscillations of the 
X-ray flux. The oscillations are aperiodic with a characteristic time scale of 14.2 days and an order of magnitude 
variation in the 2.8-20 keV flux. The oscillations may be related to optical "mini-outbursts" that have been 
observed at the ends of outbursts for other short orbital period BHC transients. The X-ray flares have durations 
between 62 and 215 s and peak fluxes that are 5-24 times higher than the persistent flux. The flares have non- 
thermal energy spectra and occur when the persistent luminosity is near 3 x 10 34 (c?/4 kpc) 2 erg s" 1 (2.8-20 keV). 
The rise time for the brightest flare demonstrates that physical models for BHC systems must be able to account 
for the situation where the X-ray flux increases by a factor of up to 24 on a time scale of seconds. We discuss the 
flares in the context of observations and theory of Galactic BHCs and compare the flares to those detected from 
Sgr A*, the super-massive black hole at the Galactic center. We also compare the flares to X-ray bursts that are 
seen in neutron star systems. While some of the flare light curves are similar to those of neutron star bursts, the 
flares have non-thermal energy spectra in contrast to the blackbody spectra exhibited in bursts. This indicates that 
X-ray bursts should not be taken as evidence that a given system contains a neutron star unless the presence of a 
blackbody component in the burst spectrum can be demonstrated. 

Subject headings: accretion, accretion disks — black hole physics: general — stars: individual 
(XTE J1650-500) — stars: black holes — X-rays: stars 



1. INTRODUCTION 

X-ray transients provide an opportunity to study accreting 
compact objects over a large range of luminosities and mass 
accretion rates. Until recently, most observations of black hole 
candidate (BHC) X-ray transients occurred when the sources 
were at very high luminosity (;> 10 36 erg s" 1 ) or very low lumi- 
nosity (<; 10 32 erg s -1 ), and much less observing time was de- 
voted to intermediate luminosities (see Chen, Shrader & Livio 
1997 for a review of observations of X-ray transients). This 
situation has improved somewhat throug h observations with 
the Rossi X-ray Timing Explorer (Bradt, Rothschild & Swank, 
1993) and also with imaging X-ray observatories. Over the 
past few years, we have been observing BHC transients dur- 
ing outburst decay with RXTE, the Chandra X-ray Observatory 
dWeisskopfet all [2002 ), and also in the radio band (see, e.g., 
Tomsick & Kaaret 2000; Corbel et al. 2001; Tomsick, Corbel 
& Kaaret 2001; Kalemci 2002). One motivation for these ob- 
servations is performing detailed studies of the transitions be- 
tween the high- soft (or s oft) and low-hard (or hard) spectral 
states (Ivan der Klis, 119951 where significant changes in the X- 
ray spectral and timing properties occur. While these tr ansitions 
typic ally occur at luminosities of about 10 36 erg s" 1 (iNowakl 
1995|, our observing program is designed to extend to much 
lower luminosities as we have recognized that new phenomena 
associated with intermediate accretion rates may be found by 
observing in this rather poorly studied luminosity regime. 

In the hard state for BHCs, it is likely that the X-ray emis- 
sion is due to inverse Comptonization of soft photons by ener- 
getic electrons, but there is uncertainty about the system geom- 
etry and the mechanism for transferring energy to the electrons. 
One possibility is that the inner edge of the accretion disk re- 
cedes from the compact object in the hard state, leaving a hot, 
quasi-spherical, optically thin region in the inner portion of the 
accretion disk. Another possibility that implies a different site 



for hard X-ray production and a different accretion geometry is 
that the hard X-ray emission is due to a large number of discrete 
regions abov e the accretion disk where magne tic reconnection 
events occur dGaleev. Rosner & Vaianalll979l) . While both of 
these possibilities describe physical processes that may occur in 
BHCs, neither provides a complete picture of accreting BHCs 
as radio observations in the ha rd state indicate the presence of a 
powerful compact jet (Fender u2001l) . The compact jet has also 
been suggested as a possible contributor to the X-ray emission 
via in verse Co mptoniz ation or possibly a synchrotron mech a- 
nism (Markof LFalcke & Fenderll2001tlMarkoff et al.Ll2003l) . 

The BHC transient XTE J 1650-500 was discovered in 2001 
September when it was detected in the X -ray band in out - 
burst by the RXTE All-Sky Monitor (ASM, Remillardl 120011) . 
Pointed RXTE observations indicated X-ray spectral and tim- 
ing properties typical of black hole systems: No X-ray pul- 
sations were detected; the energy spectrum was consistent 
with a combination of a soft component, likely from the ac- 
cretion disk, and a power-law component; and band-limited 
noise and quasi-periodic oscillations (QPOs) were observed 
(iMarkwardt. Swank & Srnithll200lHRevnivtsev & Sunvaevll2001t 
Wiinands, Mill er & Lewinll2001l). XTE J 165 0-500 was identi- 
fied in the optical (Castro-Tirado et al., 2001) and also the radio 
dGroot. Tingav. Uda lski & M illed. 1200 ll) . but the radio source 
was not resolved. XTE J 1650-500 was observed in the X-ray 
band with XMM-Newton and Chandra during its outburst. The 
XMM-Newton spectrum represents one of the best examples of 
a smeared and gravitationally red-shifted iron Ka emission fine 
that h as been detected for a stellar mass BHC dMiller et all 
2002). More recently. iHoman et alJ J2003) reported the detec- 
tion by RXTE of a 250 Hz QPO, making XTE J 1650-500 the 
sixth BHC where high frequency (>40 Hz) QPOs have been 
detected. Finally, optical observations near quiescence indicate 
a binary orbital period of 5.1 hr and an optical mass function of 
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0.64 ±0.03 M ISanchez -Fernandez et al 1 12001 . The value 
for the mass function is relatively low for a black hole sys- 
tem and could indicate that the system contains a relatively low 
mass black hole, has a low binary inclination (i < 40°), or both. 

During the 2001 X-ray outburst from XTE J 1650-500, the 
source underwent dramatic changes in its X-ray spectral and 
timing properties. Generally, the changes were not unusual for 
BHC X-ray transients. The source started off in a hard state, 
with a hard energy spectrum and a high level of timing noise, 
made a transition to a soft sp ectral state, and then made a transi- 
tion b ack to the hard state ( iKalemci et all 120031 iHoman et all 
2003|. Here, we focus on RXTE observations that were made 
well after the state transition when the source flux was between 
10 and 100 times below the transition level. In §2, we describe 
the observations and our analysis of the data. We present the 
results of the analysis in §3. The results include the discov- 
ery of X-ray oscillations and flares, and we previously rep orted 
the oscillations in an IAU Circular ( Tomsi ck et all 12002). We 
discuss the results in §4 and summarize our conclusions in §5. 

2. OBSERVATIONS AND ANALYSIS 

We obtained 97 RXTE monitoring observations of XTE J 1 650- 
500 covering the end of its 2001-2002 outburst. These "Tar- 
get of Opportunity" observations started on 2001 November 5 
(MJD 52,218) and were triggered by the gradual drop in the 
X-ray flux that was observed by the RXTE/ ASM. Of the 97 ob- 
servations, 15 occurred prior to a 28 day gap caused by point- 
ing restrictions rela ted to the source position relative to the sun. 
Kale mci et alJd2003l) report on timing analysis of these pre-gap 
observations, during which the 2.8-20 keV source rate was be- 
tween 69 and 183 s" 1 per PCU (Proportional Counter Unit), in- 
cluding counts from the top anode layer only. Here, we focus on 
the 82 post-gap RXTE observations that were made over a pe- 
riod of 182 days between 2001 December 22 and 2002 June 22. 
Immediately after the gap, the source rate had dropped to 2 1 s" 1 , 
and for the remaining observations, the source rates averaged 
over the duration of each observation are shown in Figure^. 
The RXTE observations occurred approximately every 2 days 
with the exception of the time period between MJD 52,361 and 
MJD 52,385, during which there was only 1 observation. The 
exposure times for the 82 post-gap observations were between 
384 and 3360 s with an average exposure time of 1516 s per 
observation. Long time scale (~14 day) and large-amplitude 
X-ray oscillations are apparent in the Figure^ light curve. 

For the light curve shown in Figure ^ and for much of 
the work described below, we used "Standard 2" PCA (Pro- 
portional Counter Array) data with 16 s time resolution and 
129 energy channels. We extracted light curves and energy 
spectra using scripts developed at UC San Diego and the Uni- 
versity of Tubingen that incorporate the standard software for 
RXTE data reduction (FTOOLS v5.2). We used the most re- 
cent (2002 February) release of the "Faint Source" model for 
background subtraction. Obtaining the best possible back- 
ground subtraction is important for some aspects of this work 
because of the low source count rate. To achieve this, we used 
Standard 2 data from the top anode layer only and excluded 
data from PCU 0, which has a higher background level due to 
the loss of its propane layer 5 . The light curve shown in Fig- 
ure^ indicates that we obtain high quality background sub- 
traction. For the final 8 observations, we detect little or no flux 

5 See http://lheawww.gsfc.nasa.gov/users/craigm/pca-bkg/bkg-users.html 
for a detailed analysis of the background model performance. 
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FIG. 1 .— (a) The 2.8-20 keV PCA light curve at the end of the 2001 XTE 
J1650-500 outburst. The count rate is the background subtracted, top-layer 
rate for PCU 2. The diamonds mark observations during which X-ray flares 
occurred, and the squares mark the observations where simultaneous Chandra 
data were obtained. The rates shown include the counts from the X-ray flares. 
The arrows are spaced by 14.2 days, (b) The Lomb periodogram for the PCA 
light curve. The arrow marks the most significant peak in the periodogram, 
which is at 14.2 days, but it is not a periodic signal. 

from XTE J 1650-500, and the measured flux levels are consis- 
tent with the expected leve l from the Galactic ridge emission 
dVahnia & Marshall 1 19981) in the PCA field of view for XTE 
J 1650-500. Finally, we note that for some of the results de- 
scribed below, we used a second data mode with 125 /is time 
resolution that was obtained simultaneously with the 16 s time 
resolution data. 

During 3 of the RXTE observations, we also obtained si- 
multaneous observations with Chandra (on MJD 52,298, MJD 
52,309 and MJD 52,335), and these are indicated in Figure[3i. 
Although we leave a detailed analysis of the Chandra data to a 
future paper, we refer to these observations in §3.3 and §3.5. 

3. RESULTS 

3.1. Discovery ofX-Ray Flares 

Through an inspection of the 16 s 2.8-20 keV PCA light 
curves for XTE J 1650-500, we found that bright and isolated 
X-ray flares occurred during several of the RXTE observations. 
As shown in Figure|2] some of these flares show count rate in- 
creases by more than an order of magnitude over the persistent 
flux and durations near 100 s. Even though these flares are dra- 
matic bursts of X-ray emission, significant variability is present 
in many of the other 16 s light curves. Thus, to ensure that 
we find the most extreme flares, we carried out a quantitative 
search for flares using the 16 s light curves. First, we produced 
a count rate histogram (i.e., the number of light curve time bins 
in various count rate ranges vs. count rate) for each observation. 
Then, we performed a least squares fit to these histograms with 
a Gaussian function. We obtained three parameters from the 
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fit: The Gaussian mean (r mean ); the width (a); and the ampli- 
tude. For each observation, we calculated x = (r max - r mea „)/a, 
where r mOT corresponds to the maximum count rate in the 16 s 
light curve for that observation. Based on the number of ob- 
servations and their duration, if the light curves contain only 
pure Gaussian noise, one expects approximately 1 observation 
with x > 4 and no observations with x > 5. In fact, for XTE 
J1650-500, 8 (out of 82) observations have x > 7, indicating 
the presence of a non-Gaussian component to the variability. It 
should be noted that, for light curves in general, a large value 
of x does not necessarily indicate the presence of flares because 
a long-term trend could also produce a large value of x. How- 
ever, an inspection of the light curves for the 8 XTE J 1650-500 
observations with x > 7 indicates that these light curves contain 
large and rapid increases in flux (i.e., flares). In this work, we 
focus on the 5 observations that contain the most extreme flares 
with values of x between 10 and 20. 

Closer examination of the flares in the 5 observations indi- 
cates that the properties of the flare that occurred during the 
observation on 2002 January 9 (MJD 52,283) are significantly 
different from those of the other flares. The Jan. 9 flare is ex- 
tremely hard relative to the other flares. While the peak rate 
in the 2.8-20 keV band is similar for the flares in all 5 ob- 
servations, the peak 20-60 keV count rate for the Jan. 9 flare 
is 18 ± 1 s" 1 per PCU compared to non-detections in this en- 
ergy band with upper limits of 1-2 s" 1 per PCU for the other 
flares. In fact, there is good evidence that the Jan. 9 flare did 
not originate from any source within the RXTE field-of-view 
(FOV). Based on the PCA rates observed at the peak of the 
Jan. 9 flare, a source in the FOV should have been strongly de- 
tected with a count rate of approximately 600 s" 1 in HEXTE 
(High-Energy X-ray Timing Experiment), which is the hard X- 
ray (20-200 keV) instrumen t that i s co-a ligned with the PCA 
on-board RXTE ( iRothschild et all 119981) . We examined the 
HEXTE light curves for the Jan. 9 observation, and did not find 
any evidence for a count rate increase. 

The properties of the Jan. 9 flare are very similar to 3 flares 
that were detected with RXTE during observations of other 
sources in 1998. The 1998 events have been attributed to X- 
rays and gamma-rays produced during solar flares 6 . To check 
the hypothesis that the Jan. 9 flare has a solar origin and that the 
other flares we detect during the XTE J 1650-500 observations 
do not, we obtained publicly available X-ray light curves with 
5 minute time resolution from one of the Geostationary Oper- 
ational Environmental Satellites (GOES-12), which carries the 
Solar X-ray Imager (SXI). Very bright solar flares were detected 
with SXI at the time of the 2002 Jan. 9 flare and also at the times 
of the 1998 events. However, solar flares did not occur at the 
times of the other XTE J 1650-500 flares, and we conclude that 
the Jan. 9 flare has a solar origin, but the other flares do not. 
Figure|2]displays PCA light curves with 4 s time resolution for 
the 4 XTE J 1650-500 observations (henceforth, observations 
a-d) where we detect non-solar flares. We detect 6 flares (la- 
beled 1-6 in Figure^) in these observations and report on their 
detailed properties below. 

3.2. Basic X-Ray Flare Properties 

To characterize the flares, we extracted light curves for ob- 
servations a-d with 1 s time resolution. Although PCA back- 
ground estimates are only provided with 16 s time resolution, 
we used an interpolation procedure to background subtract the 

s http://lheawww. gsfc.nasa.gov/users/keith/non_cosmic_bursts/burst_doc. html 
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FIG. 2. — The 2.8-20 keV PCA light curves for the 4 observations where 
X-ray flares occurred. The time resolution is 4 s, and the light curves are 
background subtracted. To maximize the statistics, data from all PCUs that 
were on, including PCU 0, and all anode layers were used, but the rates are 
divided by the number of active PCUs to give counts s -1 PCU~' . The flares 
are labeled 1-6, and the flare and non-flare times used for spectral analysis are 
marked with solid and dashed lines, respectively, below the light curves. The 
dates of the observations are (a) MJD 52289.7, (b) MJD 52315.7, (c) MJD 
52319.1, and (d) MJD 52358.6. 

1 s light curves. For each flare, we fitted the rise and decay with 
a function consisting of an exponential plus a constant. In most 
cases, an exponential does not provide a good description of the 
light curve if the peak of the flare is included, but acceptable fits 
are obtained when we include the portion of the light curve up 
to 85% of the flare peak. For the 12 fits (rises and decays for 
6 flares), we obtain reduced \ 2 values between 0.81 and 1.54 
with a mean reduced x 2 of 1.07, typically for 110 degrees of 
freedom. Table Ogives the e-folding rise and decay times for 
the 6 flares and also the flare durations. The durations are cal- 
culated using the exponential fits and defining the start of the 
flare to be the time when the rate rose to 20% above the con- 
stant level and the end of the flare to be the time when the rate 
reached 20% above the constant level during the decay. The 
results indicate that the flares have rise times of 3-18 s, decay 
times of 15-98 s, and durations of 62-215 s. For each flare, Ta- 
ble n alS0 includes the ratio of the peak 2.8-20 keV count rate 
in the 4 s light curves to the non-flaring (i.e., persistent) rate in 
the same band. The non-flaring light curve regions are marked 
in Figure|2] The peak rates are 5-24 times higher than the per- 
sistent rates. 

3.3. Flare and Non-Flare Energy Spectra 

We studied the spectral properties for observations a-d, by 
producing 2.8-20 keV PCA energy spectra using data from the 
flare and non-flare times shown in Figure |2] We produced en- 
ergy spectra for each flare individually and also a combined 
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spectrum with a 464 s integration time. We determined the non- 
flare times by inspection of the 16 s light curves and included 
all long, continuous segments where the count rate is relatively 
constant. We produced spectra for each of the 4 observations, 
and the exposure times for the non-flare spectra are given in Ta- 
ble |2] For spectral fitting, we used the XSPEC 11.2 software 
package. Based on fits to PCA energy spectra from observa- 
tions of the Crab nebula that were contemporaneous with the 
XTE J 1650-500 observations (see Tomsick, Corbel & Kaaret 
2001 for a more detailed description of this procedure), we in- 
cluded 0.6% and 0.3% systematic errors for energy bins below 
and above 8 keV, respectively. 

We fitted the 464 s flare spectrum with power-law and black- 
body spectral models to determine if the flare emission mech- 
anism is non-thermal or thermal. For these fits, we included 
interstellar absorption using the "phabs" model and left the col- 
umn density (A^h) as a free parameter. Using the power-law 
model results in a very good fit with x 2 / v = 23/37, while the 
result for the blackbody model is a poor fit with x 1 l v = 120/37. 
Thus, it is clear that the flare spectrum is not dominated by a 
blackbody component. Furthermore, a two-component power- 
law plus blackbody fit does not provide a significant improve- 
ment over the power-law alone, and the upper limit on the 
blackbody contribution to the two-component fit is 7% of the 
total 2.8-20 keV flux. 

The non-flare spectra are also well-described by a power-law 
model. However, for observations a and d, positive residuals 
appear between 6 and 7 keV that may indicate the presence of 
an iron Ka emission line. For observation d, where the largest 
residuals occur, the line is consistent with being narrow and 
the line energy is 6.6 ± 0.3 keV. The equivalent width (EW) is 
poorly constrained to be between 200 eV and 1.4 keV (90% 
confidence). Based on the improvement in the quality of the 
fit when the line is included in the model, an F-test indicates 
that the line is required at only the 99.4% level. For the other 
non-flare observations and for the 464 s flare spectrum the 90% 
confidence EW upper limits are a few hundred eV or greater - 
consistent with the EW for observation d. Thus, for the spectral 
fits described below, we use a model consisting of a power-law 
and a narrow iron line. 

For the flare and non-flare spectra, the power-law plus line 
spectral fits indicate that the constraints on the column density 
are poor, and we obtain 90% confidence upper limits on A^h in 
the range (2-10)x 10 22 cm" 2 . Our Chandra spectra (see §2) in- 
dicate a column density of about A^h = 6 x 10 21 cm" 2 , which 
is consistent with the PCA upper limit s, the expected inter- 
stellar value dDickev & Lockmanl 11990ft . and also the values 
inferred from X-ray and optical observations obtained during 
outburst (Miller et al. 2002; Augusteijn et al. 2001; A. Castro- 
Tirado, private communication). We refitted the spectra after 
fixing A^h to the value we obtained from Chandra, and the re- 
sults for the combined flare spectrum and the non-flare spec- 
tra are given in Table |2] While the combined flare spectrum 
does have a smaller best fit value for the power-law photon in- 
dex (r) than the non-flare spectra, it is only marginally harder 
than the non-flare spectra for observations b and d. However, 
at F = 1.59 ± 0.08 (90% confidence errors), the flare spectrum 
is significantly harder than the spectra of the two non-flare ob- 
servations with the lowest flux, a and c, for which T = 2.20^{j 26 
and r = 2.50^Q4j, respectively. This indicates that significant 
softening occurs at low flux levels, but it is not clear if this is 
related to the flaring behavior. We refitted the combined flare 



spectrum after subtracting the average non-flare spectrum. This 
results in a slightly lower (harder) value for T, but the change 
in r is less than 0. 1 . 

We also fitted the spectra for the individual flares and the val- 
ues of r obtained are given in Tabled along with the fluences 
for each flare. To determine the fluences, we first summed the 
total number of counts over the duration of the flare and sub- 
tracted the contribution from the persistent flux. We determined 
the start and end times for each flare as described in §3.2. We 
then used the spectral shape from fitting the individual flares to 
convert the total counts to a fluence in erg cm" 2 . 

3.4. Evolution ofX-Ray Flare Properties 

Figure [3] shows the light curves for the 6 flares in the hard 
(6.9-20 keV) and soft (2.8-6.9 keV) X-ray bands. While there 
are some specific light curve features that appear for only one 
of the energy bands, the main result is that there is good overall 
correlation between the soft and hard light curves. We also ex- 
amined the hardness ratio vs. time using the hard and soft light 
curves, but the only obvious trend is that the flares are some- 
what harder than the persistent emission, which is consistent 
with the spectral results reported in §3.3. For flare 6 especially, 
it is apparent that some of the specific light curve features are 
present for only one of the energy bands. During the rise of 
flare 6, near time = 400 s, there is an 8 s dip that only appears 
in the soft band light curve. Also, during the flare 6 decay (time 
= 465-540 s), oscillations with a ~20 s time scale are present 
in the hard band light curve only. In addition, flares 2 and 3 
show some evidence for variability in the hard band that is not 
present in the soft band. 

3.5. Long Time Scale Oscillations Observed with RXTE and 

Chandra 

Here, we focus on the long time scale and large amplitude 
X-ray oscillations seen in Figure^. To determine the charac- 
teristic frequencies of the oscillation and to search for period- 
icities, we produced a Lomb periodogram ( Pre ss et a u cm 
of the time series, and it is shown in Figure^). The highest 
amplitude peak in the periodogram is at a frequency of 0.070 
days" 1 (14.2 day period), but it is about a factor of 2 below the 
3-ct detection limit for a periodic signal. The fact that there is 
excess power from 0.06 to 0.09 days" 1 suggests the presence 
of an aperiodic signal, and inspection of the light curve agrees 
with this conclusion. In Figure [^, the arrows are spaced by 
14.2 days, and while the first 4 peaks are consistent with this 
period, the agreement does not continue after MJD 52,330. 

As the PCA is not an imaging instrument and has a rela- 
tively large field-of-view (1° radius), it is necessary to consider 
whether another source besides XTE J 1650-500 could be re- 
sponsible for the phenomena detected by the PCA. Thus, it is 
important that the Chandra observations mentioned in §2 firmly 
establish that XTE J1650-500 is the source of the oscillations. 
For all three Chandra observations, a source is strongly de- 
tected at the position of the X TE J 1650-500 radio counterpart 
dGroot. Tingav. Udalski & Millerl 1200 lft . The Chandra spec- 
tra are well-described by a power-law with interstellar absorp- 
tion, and by fitting the spectra, we derive 0.5-10 keV fluxes 
of 5 x 10~ n erg cm" 2 s" 1 for the first two observations and 
6 x 10~ 12 erg cm" 2 s" 1 for the third observation. Comparing 
these fluxes to the PCA count rates shown in Figure^ indi- 
cates that the flux levels measured by Chandra and the PCA are 
well-correlated. As the Chandra observations prove that XTE 
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FIG. 3. — Light curves in 2 energy bands for the 6 X-ray flares. The time 
resolution is 4 s, and the light curves are background subtracted. To maximize 
the statistics, data from all PCUs that were on, including PCU 0, and all anode 
layers were used, but the rates are divided by the number of active PCUs to 
give counts s PCLT 1 . For each X-ray flare, the soft 2.8-6.9 keV band light 
curve is shown on the bottom. The hard 6.9-20 keV band light curve is shown 
on top and has been offset by 10 counts s for clarity. 



J 1650-500 was still active when the flares were observed, it is 
very unlikely that the flares come from another source. Still, we 
used the SIMBAD database to search for X-ray sources in the 
PC A FOV, which is a 1° circle centered on the XTE J 1650-500 
position (I = 336.70°, b = -3.44°). We did not find any sources, 
such as X-ray binaries, that are likely to contribute significantly 
to the X-ray flux detected by the PCA. Further, the non-thermal 
flare spectrum rules out the possibility that the flares come from 
an X-ray burster in the PCA FOV (see §4.5 for more discussion 
of X-ray bursts). 

4. DISCUSSION 

During RXTE observations of the BHC XTE J 1650-500, we 
have discovered 6 X-ray flares lasting from 62-215 s and peak- 
ing at count rate levels between 5 and 24 times the persistent 
rates. The flares occurred at the end of the first and currently 
only X-ray outburst that has been detected from this X-ray tran- 
sient when the persistent flux was a factor of 10 3 below the the 
peak flux level for the outburst. Although strong X-ray vari- 
ability is commonly observed for black hole X-ray binaries, the 
isolated X-ray flares that we observe are unusual. In addition, 
we find long term aperiodic oscillations with a characteristic 
time scale of 14.2 days. Below, we discuss the implications of 
the XTE J 1650-500 oscillations and X-ray flares by compar- 
ing them to phenomena that have been seen for other sources, 
including long and short time scale variability in black hole bi- 
naries, the X-ray flares that have been detected from Sgr A*, 
the super-massive black hole at the center of the Galaxy, and 
X-ray bursts in neutron star binaries. 



4.1. XTE J 1650-500 Oscillations 

The fact that the ^14 day X-ray oscillations are an aperi- 
odic phenomenon likely rules out the possibility that they are 
directly related to disk precession, and they are obviously not 
orbital since it is now known that XTE J 16 50-500 has an or- 
bital period of 5.1 hr ( Sanchez-Ferna ndez et all f2002l) . How- 
ever, oscillations have previously been observed in the optical 
at the ends of outbursts for a number of short orbital p eriod 
BHC transients as well as dwarf novae ( lKuulkersl[T998l) . The 
general properties of the optical oscillations can be explained in 
the framework of disk instability m odels (DIM) f or dwarf novae 
dOsaki. Shimizu & TsugawaLll997tfHameurv . Lasota & Warneij 
2000) and possibly for X-ray transients dMenouet allT loOO). 
For DIM, the level of emission increases when the disk be- 
comes hot enough for the hydrogen in the disk to be ion- 
ized. This leads to a higher viscosity and an increase in the 
mass accretion rate onto the compact object. Oscillations can 
be produced by the propag ation of heating and then cooling 
waves through the disk (Hameury. Lasota & Warner, 2000). 
For BHC systems, the best examples of these optical oscil- 
latio ns (also called "mini-outbursts") include GRO J04 22+32 
dCastro-Tirado et allll993[IChevalier & Ilovaiskvl[l995l) . GRS 



1009-45 dBailvn & OrosTll995h . and, more recently. XTE 
J1859+226 lZurita et all 120021) . For XTE J 1650-500, optical 
oscillations occurred contemporaneously with the X-ray oscil- 
lations we are reporting in this work (C. Bailyn, private com- 
munication). The optical oscillations have similar characteristic 
time scales of 10-100 days for all 4 sources. While X-ray cover- 
age during optical oscillations was not obtained for GRS 1009- 
45 and was sparse for GRO J0422+32 and XTE J 1859+226, 
the X-ray observations obtained for the latter 2 sources showed 
that they were active in the X-ray band during the oscillations 
dShrader et all Il997t iTomsick & Heindll 120001) . Another short 
orbital period BHC transient, 4U 1543-47, showed large varia- 
tions in the X-ray band in a poorly sampled light curve obtained 
at the end of its 1971-1972 outburst dChen. Shrader & Liviol 
119971) . Overall, the similarities between the XTE J1650-500 
oscillations and those that have been seen for other short orbital 
period BHC transients suggests that they are related phenom- 
ena. If this is the case, our observations of XTE J 1650-500 
represent the first time that good X-ray coverage has been ob- 
tained during a series of mini-outbursts. 

4.2. XTE J 1 650-500 X-Ray Flares and BHC Variability 

While it is likely that the XTE J 1650-500 oscillations are 
similar to phenomena that have been observed before (mostly 
in the optical band) for BHCs, it is not clear that the same can 
be said for the X-ray flares. In certain spectral states, BHC 
light curves certainly do exhibit a high level of variability, and 
their X-ray lig ht curves conta in flares that have inspired "shot 
noise" models <Terrelllfl972l However, in these models, it is 
the combination of a very large number of shots that can lead to 
the overall variability observed in BHCs, which can be as fast 
as millisecond time scales. The shots, and also the flares that 
are typically observed in BHC light curves, have much smaller 
amplitudes relative to the average flux levels and much shorter 
time scales than the XTE J1650-500 flares (see Miyamoto et 
al. 1992 for examples of light curves for other BHCs). 

We found one example in the literature of a black hole sys- 
tem that exhibited a series of X-ray flares with some similarities 
to the XTE J1650-500 flares. V4641 Sgr is an unusual X-ray 
transient that was discovered in 1999 February by RXTE and 
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BeppoSAX. It is often referred to as a "fast" transient because it 
had an outburst in 1999 September that only lasted for about a 
day rather than the typical several months. The source showed 
extreme X-ray variability when it was bright, and, at the end of 
its most active period in 1999 Septem ber, it exhi bited relatively 
isolated X-ray flares ( Wiinands & van derKlisi 12000). Three 
or four X-ray flar es occurred in a 200 s lig ht curve shown in 
Figure lc of Wiina nds & van der Klisl (120001) . The X-ray flares 
had rise times <10 s, durations of 20-50 s, flux increases by 
factors of 10-50 over the non-flaring flux, and roughly expo- 
nential decays at the ends of the flares. The similarities of the 
flare parameters and the fact that the flares for both V4641 Sgr 
and XTE J 1650-500 occurred when the non-flaring flux was 
several orders of magnitude below the peak outburst flux for 
the respective sources leads us to suggest the possibility that 
the flares from the two sources could be related. 

Although we do not believe that similar X-ray flares have 
been seen before in BHC sources other than XTE J 1650-500 
and possibly V4641 Sgr, this could partly be a selection ef- 
fect. As discussed in §1, BHCs have been most often ob- 
served in outburst at high luminosity or in quiescence. From the 
optical observations of XTE J 1650-500 that have been made 
dSanchez-Fernandez et alll2002tlAugusteiin. Coe & Grooil200ll 
iGarcia & Wilkesl 1200 21 we estimate that the source distance is 
between 2 and 6 kpc, and here and below we use a fiducial 
value of 4 kpc for luminosity estimates. Thus, the typical non- 
flare X-ray flux of 1.7 x 10" 11 erg cm" 2 s" 1 (2.8-20 keV, un- 
absorbed) corresponds to an luminosity of 3 x 10 34 (c//4 kpc) 2 
erg s" 1 , which is in a luminosity regime where BHCs have not 
been as well-studied. Thus, we cannot rule out the possibil- 
ity that X-ray flares commonly occur for BHCs at intermediate 
luminosities. 

The fact that multiple X-ray flares are detected from XTE 
J 1650-500 during observations a and b, suggests the possibility 
that XTE J 1650-500 entered into a state (or perhaps an accre- 
tion regime) that favors flare production. It is interesting that 
the flares for observations a and b are both separated by 800- 
900 s, and this time scale may provide a clue about the flare 
emission mechanism. The observation d light curve also sup- 
ports the idea that XTE J 1650-500 enters into a flaring state 
since a short flare is seen at the end of that observation (near 
time = 1850 s). Furthermore, the RXTE observation follow- 
ing observation d (on MJD 52361) is one of the observations 
for which higher variability was detected, and short flares are 
also detected during that observation. Finally, a higher level of 
variability was detected for the RXTE observation that occurred 
between observations b and c on MJD 52318. 

4.3. X-Ray Flares and Black Hole Accretion Theory 

The XTE J 1650-500 X-ray flares demonstrate that physi- 
cal models for BHC systems must be able to account for the 
situation where the X-ray flux increases by a factor of up to 
24 on a time scale on a time scale of seconds (based on flare 
5). While this time scale does not provide a useful model- 
independent constraint on the size of the emission region be- 
cause the physical size of the system (i.e., the binary separa- 
tion) is only a few light-seconds, we have considered whether 
the flares have implications for certain accretion scenarios. For 
flares produced by a change in mass accretion rate, the flux is 
expected to change on the viscous time scale in the X-ray emit- 
ting region of the accretion disk. For thin (H <C R) and thick 
(H ~ R) accretion disks, where H is the accretion disk thick- 
ness and R is the radial distance from the compact object, the 



viscous time scale is given by t v = a (R/H) 2 ta, where a is 
a dimensionless parameter that is proportional to the viscosity 
and td is the dynamical time scale ( Frank. King & Rainelfl992l 
The dynamical (i.e., Keplerian) time scale depends on the mass 
of the compact object, M, and R = rRs (Rs = 2GM/c 2 , where G 
is the gravitational constant and c is the speed of light) accord- 
ing to td = 88/is (M/M©) r 3//2 . Based on the current distribution 
of black hole masses, it is likely that XTE J 1650-500 contains 
a black hole with a mass of 5-15M Q , a nd we assume a mass 
of_10M Q . In the extreme case of a = 1 (tFrank. King & RaineL 
fl992l) and a thick accretion disk (R = H), f,, = 1 s corresponds 
to r = 1 10, while a more realistic value of a = 0.1 gives r = 23. 
Thus, if the flares are due to an increase in the mass accre- 
tion rate, they must originate in a region smaller than the inner 
^20-100 Rs of the accretion disk. While this does not strongly 
constrain theoretical models, the possibility of larger hard X- 
ray emission regions has been previously suggested (e.g., Esin, 
McClintock & Narayan 1997). 

Some theoretical models for accreting black hole systems 
predict highly variable X-ray emission at low luminosities. 
iMerloni & Fabianl (120021) discuss the possibility of a "coronal 
outflow dominated accretion disk model," where a geometri- 
cally thin and optically thick accretion disk dissipates much of 
its gravitational energy in a magnetic corona. This is an attrac- 
tive model since it incorporates angular m omentum transport 
via magnetic fields ( Balb us & HawlevtfT99H and also the pres- 
ence of outflows such as the compact jets that have been ob- 
serve d from accreting black holes at low luminosities dFendeil 
2001). In this model, non-thermal X-ray flares can be pro- 
duced through gradual magnetic field gene ration that c an lead 
to a brightening of the corona ( Merloni & Fabian, 2002). In an- 
other model that has been suggested for low luminosity black 
hole systems, the black hole accretes via a "hot settling flow" 
dMedvedev & Murravtl2002T) . Here, the black hole is magnet- 
ically coupled to the accretion flow, allowing a rapidly rotat- 
ing black hole to trans fer its angular mom entum to the accre- 
tion flow as in Blandford & Znaiek ( 1977). Hard X-ray spec- 
tra and high varia bility are expected for this accretion mode. 
In fact, Medvedev & Murrav (2002) specifically predict that 
XTE J 1650-500 may exhibit this behavior based on the spec- 
tral d etermination that it contains a rapidly rotating black hole 
dMilleret alll2002l) . 

4.4. Comparison between the X-Ray Flares in XTE 
J1650-500andSgrA* 

Here, we consider the possibility that the XTE J 1650-500 
flares are related to the isolated X-ray flares that have been 
detec ted from Sgr A* (Bag anoff et all 120011 iGoldwurm et all 
2003). Such a comparison is interesting because these are both 
low luminosity BHC or black hole systems, and the overall 
shapes of the flare light curves and the evolution of the energy 
spectra are similar in the two cases. For the Sgr A* flare seen by 
Chandra in 2000 October, the flux increases by a factor of about 
45, and during t he 10 ks flare, the X-ra y spectrum hardens and 
is non-thermal (Baganoff et al., 2001). While the results for 
the 2000 October flare have been reported in the most detail, 
subsequent Chandra observations of Sgr A* yielded the detec- 
tion of several more flares, and these suggest that flux increases 
by a factor of 5-25, si milar to the XTE J 16 50-500 flares, are 
probably more typical (Baganoff et all 120021) . While the flare 
durations and rise times are about a factor of 100 less for XTE 
J 1650-500 (durations of 100 s for XTE J 1650-500 compared 
to 10 ks for Sgr A* and rise times of 10 s for XTE J 1650-500 
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compared to 1 ks for Sgr A*), it is notable that that the ratios 
for the two quantities are approximately the same. For Chan- 
dra observations of Sgr A*, the power-law photon index (r) is 
2.2^Q7 duri ng non-flare times an d 1.3^g during the 2000 Oc- 
tober flare (Baganoff et III 12.001), which is consistent with the 
spectral evolution we observe for XTE J 1650-500. 

Although the phenomenology of the flares in Sgr A* and 
XTE J 1650-500 have similarities, we must also consider the 
differences between these systems. Beyond the fact that the 
mass of the Sgr A* black hole is 2.6 x 10 6 M Q dEckart & Genzell 
119971) compared to the likely ~10M Q black hole in XTE 
J 1650-500, the flares we observe for XTE J 1650-500 oc- 
cur when the persistent Eddington-scaled luminosity is much 
higher than for Sgr A*. For a distance of 4 kpc and a compact 
object mass of 10M Q , the persistent flux of 10~ n erg cm" 2 s" 1 
for XTE J1650-500 (see TableQ gives LjL Edd = 2 x 10" 5 . For 
Sgr A*, the persistent X-ray luminosity is 2.2 x 10 33 erg s -1 
( Baganof f et all 120011) . indicating L x /L E dd = 5 x 10" 12 . Even 
if we conside r the bolometri c luminosity of ~40 37 erg s" 1 for 
Sgr A* dNaravan etall l 19981). L/L Edd = 2 x 10" 8 , which is three 
orders of magnitude lower than XTE J 1 650-500. Thus, we con- 
sider the possibility that the flares in the two systems are related 
with the caveat that their accretion properties are likely to be 
significantly different. 

If the flares are produced by similar processes, the fact that 
the flare time scales are longer for Sgr A* than for XTE J 1650- 
500 is not surprising given the much larger black hole mass for 
Sgr A*. However, it is not obvious why the time scales would 
be different by only a factor of 100. For example, dynami- 
cal and viscous time scales at the same accretion disk radius, 
r = R/Rs, are proportional to the black hole mass. Thus, if the 
flares were produced by an increase in th e mass accretion rate , 
the viscous time scale would be relevant (Liu & Melia, 2002), 
and a diffe rence in time scale by a factor of 3 x 10 5 would be 
predicted. Bagan off et alJ ( 1200 ll) discuss the possibility that 
the Sgr A* flare emission could have a synchrotron or a syn- 
chrotron self-Compton (SSC) mechanism. This could occur if 
electrons were accelerated in magnetically active regions of an 
accretion flow (perhaps due to magnetic reconnection events). 
In these cases, the flare durations may be related to the syn- 
chrotron or SSC cooling time scales, which a re proportional 
to 7 -1 ZT 2 (Begelman, Blandf ord & Ree4 fl984). where 7 is the 
Lorentz factor for the X-ray producing electrons and B is the 
magnetic field strength. The fact that the flare X-ray spectra 
are similar for Sgr A* and XTE J 1650-500 may suggest that 
7 is not extremely different for the two systems, but, if B was 
an order of magnitude larger in the XTE J 1650-500 emission 
region, then the predicted time scale would be a factor of 100 
lower for this system, as obs e rved. However, the calculations of 
Idi Matteo. Celotti & Fabianl ( 1 19971) suggest that the difference 
in magnetic field strength should be much greater than an order 
of magnitude. Finally, jet-based models have be en suggested to 
explain the Sgr A* flares (Marko ff et alll20 01). In these mod- 
els, the X-ray emission has a synchrotron or SSC mechanism, 
but the electrons are continuously shock accelerated so that the 
cooling time scales do not relate to the predicted flare duration 
in a simple way. More work is required to determine if any of 
these models can explain the time scale ratio of 100 between 
the XTE J1650-500 and Sgr A* flares. 

4.5. X-Ray Flares and X-Ray Bursts 

Accreting low magnetic field neutron stars exhibit type I X- 
ray bursts due to thermonuclear instabilities that result from ac- 



cretion onto the neutron star surface ( iLewin. van Paradiis & TaamL 
fl995l) . While XTE J1650-500 is a BHC and is not expected to 
produce type I X-ray bursts, the flare rise times, the decay times, 
and the ratio between the peak X-ra y flux and the persistent flux 
are al l typical of type I X-ray bursts ( ILewin. van Paradiis & TaamL 
1995). Even though the light curve shapes for the flares are typ- 
ical of type I X-ray bursts, the spectral evolution is not. The 
emission mechanism for type I X-ray bursts is thermal, and X- 
ray bursts have energy spectra that are dominated by a black- 
body component. As shown in §3.3, the possibility that the 
XTE J 1650-500 X-ray flares have a significant blackbody com- 
ponent is strongly ruled out. 

Although type II X-ray bursts are not as well understood and 
have only been detected from two sources to date, it is notable 
that some of the XTE J 1650-500 flare light curves have prop- 
erties similar to type II bursts. For flare 6, the variability lead- 
ing up to the flare (time = 0-300 s in Figure |2ji, the dip right 
before the flare (time near 360 s), and the oscillations during 
the decay have similarities to the light curves of type I I bursts 
(Lewin, van Paradiis & TaamL fl995tlLewin et allfl996l) . How- 
eve ^ like type I bursts, type II bursts als o exhibit therma l spec- 
tra ( ILewin. van Paradiis & Taam, Il995t iMarshall et all [2001) 
in contrast to the XTE J 1650-500 flares. If XTE J 1650-500 
contains a black hole, then our conclusion that the X-ray flares 
are not type I or type II X-ray bursts indicates that X-ray bursts 
should not be taken as evidence for a neutron star in the system 
unless it can be shown that the burst energy spectrum is dom- 
inated by a blackbody component. This conclusion may have 
implications for estimates of the numbers of neutron star bina- 
ries in the Galaxy if a fraction of the putative X-ray bursters 
without high quality burst spectra actually contain black holes. 

While the energy spectrum provides strong evidence against 
the XTE J 1650-500 X-ray flares being X-ray bursts, we note 
that this possibility is also disfavored from luminosity consid- 
erations. Typically, peak X-ray burst luminosities are between 
0.1 L E dd an d L E dd- However, the flux at the peak of the XTE 
J 1650-500 outburst is about 100 times higher than the peak of 
the strongest X-ray flare, indicating that if the X-ray flares are 
close to L E dd, then the source was super-Eddington for most 
of the outburst, which seems implausible. In addition, given 
the observed flux, a peak X-ray flare luminosity of 0. \L E dd im- 
plies a distance between 28 kpc and 89 kpc for compact ob- 
ject masses between 1 M and 10M Q , and these large distances 
are ruled-out by the optical detections of the source near quies- 
cence. However, these luminosity arguments may not be abso- 
lute proof against X-ray bursts because there are some unusual 
cases where bursts in other systems have been classified as type 
I X-ray bursts even though they are significantly sub-Eddington 
(e.g., Gotthelf & Kulkarni 1997; Garcia & Grindlay 1987). On 
the other hand, for the examples in the two preceding refer- 
ences, it has not been clearly demonstrated that non-thermal 
burst energy spectra are ruled out. 

5. SUMMARY AND CONCLUSIONS 

In this work, we report on RXTE observations of the BHC 
transient XTE J 1650-500 made at the end of its first, and cur- 
rently only, outburst. By observing at low luminosities, we 
found new or previously poorly studied X-ray phenomena, in- 
cluding bright ~100 s X-ray flares and long time scale oscilla- 
tions of the X-ray flux. Over several months, the X-ray light 
curve shows strong aperiodic oscillations with a characteris- 
tic time scale of 14.2 days during which the X-ray flux varies 
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by an order of magnitude. We argue that the oscillations are 
probably related to the optical oscillations (also called "mini- 
outbursts") that have been observed for other short orbital pe- 
riod BHC transients such as GRO J0422+32, GRS 1009-45, 
and XTE J 1859+226. If this is the case, our observations repre- 
sent the first time that good X-ray coverage has been obtained 
during a series of mini-outbursts. 

A uniform analysis of the 16 s RXTE light curves leads to 
the detection of 6 X-ray flares in 4 observations with dura- 
tions between 62 and 215 s and peak fluxes that are 5-24 times 
higher than the persistent flux. The flare energy spectra are non- 
thermal, being well-described by a power-law with a photon in- 
dex of r = 1 .59 ±0.08 (90% confidence errors). The upper limit 
on a blackbody contribution to the flare is 7% of the 2.8-20 keV 
flux. The flares are isolated events that occur when the persis- 
tent luminosity is low, near 3 x 10 34 (c//4 kpc) 2 erg s" 1 (2.8- 
20 keV), making them unusual, if not unique, for BHCs. The 
XTE J 1650-500 X-ray flares demonstrate that physical mod- 
els for BHC systems must be able to account for the situation 
where the X-ray flux increases by a factor of up to 24 on a 
time scale of seconds (based on flare 5). We briefly discuss two 
theoretical models for accreting black holes that predict h ighly 
variable X-ray emission at low lum inosities (Merloni & Fabian, 
l2002tlMedvedev & Murravll2002l) . 

We also compare the XTE J 1650-500 flares to the X-ray 
flares that have recently been observed from Sgr A* by Chan- 
dra and XMM-Newton as these are both accreting black holes 
or BHCs that exhibited isolated X-ray flares during periods of 
low persistent luminosity (although we note that L x /Lgdd and 
bolometric L/LEdd are considerably lower for Sgr A*). The 
increases in flux and spectral evolution for the flares in the two 
systems are comparable, while the time scales (rise time and du- 
ration) are about 100 times faster for XTE J 1650-500. Given 
the factor of ^3 x 10 5 difference in compact object mass be- 
tween the two systems (assuming that XTE J 1650-500 contains 
a 10M Q black hole), it is not surprising that the time scales 
are faster for XTE J 1650-500. However, the simplest physical 
models suggest that the time scale difference should be consid- 
erably larger. Overall, we cannot rule out the possibility that the 
flares in the two systems are related, and more work is required 
on this topic. 

Finally, we note that some of the XTE J1650-500 flares have 
light curves with similarities to type I X-ray bursts that are seen 
in neutron star systems, including rise times of seconds and ex- 
ponential decays with e-folding times of tens of seconds. How- 
ever, the non-thermal energy spectra for the flares are in contrast 
to the spectra of type I (and also type II) X-ray bursts, which are 
dominated by blackbody emission. Thus, the XTE J 1650-500 
flares indicate that X-ray bursts should not be taken as evidence 
for the presence of a neutron star in a given system unless the 
presence of a blackbody component can be demonstrated. 
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Table 1 
Basic Flare Properties 



Flare 


MJD n -52,000 Rise Time* 


Decay Time* 


Duration* Count Rate 






Number 


(days) (s) 


(s) 


(s) Increase c 




Fluence e 



1 289.69046 7±3 25±6 65 ± 1 1 12.4±1.4 1 ^ 2 -Q2j 2.8xl(T 

2 289.70009 4.2 ±1.6 17±4 62±7 8.7± 1.3 l-62+|| 2.7 x 10" 9 

3 315.65935 13 ±4 15 ±4 67 ±9 8.2 ±1.1 l-5l3 : i 2.5x10-' 

4 315.67009 18 ±10 33±11 118±24 5.8±1.0 1.63^ : |= 2.9 x lO" 9 

5 319.11676 2.9±0.6 37±5 75±8 23.5± 1.8 1.68^o : |i 3.9 x lO" 9 

6 358.59491 18 ± 3 98±8 215 ±14 15.6±1.3 1.61^;|| 1.2xKr 8 

"Modified Julian Date at the peak of the X-ray flare. 
"See text for precise definition. 

r Ratio of the peak count rate in the 4 s light curves to the persistent 2.8- 
20 keV rate. 

rf Power-law photon index for the flare energy spectra with 90% confidence 
errors. 

c 2.8-20 keV burst fluence in erg cm~ 2 for the duration of the flare. 



Table 2 

Flare and Non-Flare Energy Spectra 



Obs. 


F/NF" 


Exposure" (s) 


r c 


Flux" 


a-d 


F 


464 


1.59 ±0.08 


1.09 x 10" 1U 


a 
b 


NF 
NF 


1392 
1712 


2.20 +031 


1.47 x 10" 11 
1.67 x 10-" 


c 
d 


NF 
NF 


1056 
1088 


2 cn+ftS 
z - 3u -g.43 


9.61 x 10~ 12 
1.69 x 10-" 



"Flare (F) or non-flare (NF). 
"Exposure time included in the spectra. 

c Power-law photon index for the energy spectra with 90% confidence errors. 
d 2.8-20 keV absorbed flux in erg cirr 2 s _1 . 



